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There  is  a  mistake  in  Farhat’s  model  published  in  volume 
138 

There  is  a  mistake  in  Farhat’s  model  and  results  [1].  The 
primary  goal  of  this  letter  is  to  disprove  his  model  and  results. 
In  his  model,  wrong  expression  of  Eqs.  (3)  and  (4)  induced  the 
mistake  of  the  modified  Butler- Volmer  equation.  According  to 
the  article,  kf  and  kb  are  rate  constants  (cm-1  s-1),  Cr  and  Cp 
are  concentrations  (mol  cm-3),  S  is  active  surface  area  (cm2). 


Then  the  units  of  k(CRS  and  kb CpS  are  mol  cm 


-4  „-l 


,  but  the 


units  of  Vf  and  nj,  are  m°l  cm-2  s- 1 .  Therefore,  Eqs.  (3)  and  (4) 
[1]  is  wrong.  Now,  we  give  the  corrected  equations  and  deduc¬ 
tion,  and  definitions  of  symbol  are  held  as  the  same  of  Farhat’s 
model. 

According  to  the  half  cell  electrode  reaction  Eq.  (2)  [1],  the 
forward  and  backward  reactions  rates,  respectively,  are  as  fol¬ 
lows: 


Vf  =  kfCR  = 


itb  =  kbCp  = 


If 

rrFS 

h 

nFS 


(1) 


(2) 


The  net  reaction  rate  is,  therefore,  the  difference  between 
forward  and  backward  reaction  rates,  i.e.: 


Pnet  —  Vf  —  Vb  =  kfCR  -  kb  Cp  = 


If  -  h 
nFS 


rearranging: 

I net  =  If  ~  h  =  nFS(k[CR  -  kbCp) 


(3) 


(4) 


According  to  Eqs.  (7)  and  (8)  of  Farhat’s  model  [1],  we  get: 


where: 
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As  can  be  seen,  Eq.  (8)  is  quite  different  from  Eq.  (13)  in 
Farhat’s  article.  According  to  Eqs.  ( 14)— ( 17)  in  Section  3.1  in 
Farhat’s  article,  we  get: 
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where  mJA  is  the  catalyst  loading  (g  cm-2),  which  is  a  important 
parameter  of  PEMFC  catalyst  layer.  The  main  control  equation 
should  be  changed  into: 
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In  conclusion,  the  influence  of  different  microstructural 
parameters  in  a  catalyst  layer  should  be  studied  by  examining 
model  predictions  given  by  Eq.  (10)  but  not  Eq  (24)  of  Farhat’s 
model. 
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where: 


7o  =  nFSk0C^  “}Cp" 


(6) 


Dividing  by  the  nominal  area  A  and  rearranging,  then  the 
activation  polarization  for  the  oxygen  reduction  reaction  at  the 
cathode  will  become: 
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